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Abstract. The chaotic nature of a storage-ring Free Electron Laser (FEL) is investigated. The derivation 
of a low embedding dimension for the dynamics allows the low-dimensionality of this complex system to 
be observed, whereas its unpredictability is demonstrated, in some ranges of parameters, by a positive 
Lyapounov exponent. The route to chaos is then explored by tuning a single control parameter, and a 
period-doubling cascade is evidenced, as well as intermittence. 

PACS. 05.45.-a Nonlinear dynamics and nonlinear dynamical systems - 82.40.Bj Oscillations, chaos, 
and bifurcations - 05.45.Pq Numerical simulations of chaotic systems - 42.65.Sf Dynamics of nonlinear 
optical systems; optical instabilities, optical chaos and complexity, and optical spatio-temporal dynamics 
- 41.60.Cr Free-electron lasers - 29.20.Dh Storage rings 

1 Introduction chaos frequently occurs in these systems, and links have 

been done between FEL chaos and plasma chaos jSiiSj. 

Perturbations of systems have been shown to generically 

lead to bifurcations in the dynamics, for example in biol- in the case of FELs, a thorough understanding of the 
ogy [T], ecology [5], chemistry [3] or astrophysics |4j. The nonlinear dynamics of the system is of paramount impor- 
bifurcations is eventually replaced by chaos, which has tance, since their ultimate goal is to deliver stable light 
been observed in a wide range of plasma/wave interaction pulses for e.g. chemistry, biology and surface studies ex- 
such as plasma systems [S] and FELs t6iiTj. The route to periments [TUllTTl fT^. In particular, the Sensitivity to Ini- 
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can drive the system into very different dynamical regimes, model. It was shown that the recorded temporal series 
even for tiny variations of the initial state. Storage Ring were chaotic and assumption was made that, as for con- 
FEL (SRFEL) are particularly rich from a nonlinear dy- ventional lasers, the route to chaos goes through a period 
namics point of view. The evolution of the pulse train doubling cascade [23j. Afterward, a more complete model, 
shape from stationary states to limit cycles 13 has been which integrated the presence of an asynchronism between 
demonstrated. The transition from a stationary state to light and electron pulses (called detuning), was proposed 
a limit cycle arises through a Hopf bifurcation [Ml . Using [M] and employed to study the presence of chaos in the 
control of chaos methods, these limit cycle regimes were dynamics. It revealed that a period-doubling bifurcation 
suppressed by stabilizing the coexisting unstable station- can arise under changes in two control parameters, 
ary state [TSlfHlfTCj . For example, the ELETTRA FEL We here characterize this chaotic dynamics by deriving 
needs this type of control because of the lack of natural the embedding dimension and the Lyapunov exponent of 
stationary state. A topological approach has also being the dynamics. Then, a single control parameter is finely 
investigated and stable orbits have been identified [17]. tuned to explore the route to chaos, and to reveal the 
FEL nonlinear studies have also enabled to highlight in- period-doubling cascade, as well as the presence of inter- 
tensity holes in the spectrum associated to drift insta- mittence. 
bilities |18J in the spectro-temporal diagram of the FEL 

intensity distribution. More generally, a SRFEL is a sys- 2 FEL temporal dynamics 

tem which presents a localized structure in presence of ad- 

During the propagation through a periodic permanent 

vection. As a consequence, this phenomenon observed in 

magnetic field structure (an undulator) , a relativistic elec- 

SRFEL systems is expected to affect a range of systems 

tron beam emits an electromagnetic wave at its funda- 

displaying localized solutions with either local or global 

mental wavelength and its odd harmonics [25j . The light 

coupling. The drift instabilities can be stabilised thanks 

amplification results from the interaction between the rel- 

to an optical feedback inspired from the community of 

ativistic electron beam in the undulator and the store opti- 

conventional laser [TOll^H] , a technique already tested on 

cal pulse in the optical cavity. The laser gain [55] depends 

the UVSOR FEL YT]. 

both on the electron beam and the undulator character- 
Here, we are particularly interested in the chaotic fea- istics. Amplification takes place at a wavelength close to 
tures of the dynamics and their conditions of emergence, the fundamental, tunable by a modification of the mag- 
Chaos studies have first been performed on the Super- netic field of the undulator. The gain reaches the satura- 
ACO and AGO storage ring FEL thanks to an external tion level defined by cavity total losses at the expense of 
modulation [21 and explained with a three-dimensional an increase of the energy spread of the electron bunch. 
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The laser is pulsed at the repetition frequency of the relaxation dynamics. The spontaneous emission ig repre- 

electron beam in the undulator (in the MHz range), and sents the initial noise from which the laser starts. The 

the length of the optical resonator must be carefully ad- equilibrium solution of the equations (|lap and (|lb[) cor- 

justed to satisfy the longitudinal synchronization between responds to / = 1 and E = 1. For a small perturbation 

the optical pulses and the electron bunches. In practice, on around / and S and by neglecting the noise, the equation 

storage rings FEL, the synchronization condition is con- 1 can be linearised, leading to a pulsed laser train shape 

trolled by the frequency of the Radio-Frequency (RF) cav- at a resonant frequency: 



ity, since it modifies the orbit length of the electrons. The p _ 1 Gq — L 

^' -ttV 2Tor, ■ 

temporal detuning S is the difference between the round 

The model given by equations (jlap and (jlbp repro- 

trip period of the light pulse in the optical cavity and 

duces well the presence of a constant and a pulsed pulse 

the electron bunch spacing Tq. A perfect synchronism is 

train shape of the laser. Yet it considers only the pulse 

required to obtain a stable stationary state of the laser 

train shape, and does not take into account the tempo- 

(pulsed at the MHz range) [18], for which the pulse train 

ral mismatch between the light and electron pulses, the 

shape is constant. Nevertheless, beyond a given threshold 

so-called detuning S. To modelize this detuning, a longi- 

of detuning, the laser exhibits a pulsed shape in the mil- 

tudinal distribution for the intensity has to be introduced, 

lisecond range. This can be simply understood with the 

which describes the pass-to-pass interaction process 
model describing the coupled evolution of the laser pulse ^ 

, y„(r) = (l-L)[l + G„(T)]y„_i(T)+z,expf-^)3;a) 

intensity 1 (pulse tram shape) and oi the normalized en- \ ^c^^ / 

ergy spread E of the electron bunch [57]: ^„ = ^n-i [^n-i — / yn-i{T)dT], (3b) 

''"s J -oo 

fiT J where w„ stands for the longitudinal distribution (centered 

-^-[G,~L){l-S) + ^, (la) 

, „ „ on the center of mass of the electron bunch) of the laser 

^ = (lb) . ... 

'''' '''s pulse at the n-th pass in the gain medium (electron bunch 

.,, ^ / 2 2\ 1 1 2 2\ 1 • ii 1 in the undulator). This bunch has a normalized energy 

with L = [a — crQ)/(crg — CTq), a being the energy spread, ' 

J i ii 1 i i spread S at the 7i-th pass and a rms longitudinal width 

(To the energy spread at the laser start up, CTe the energy ^ o 

1 , , , , ; .ru . i- i- rr ■ ,1, cFt, supposed to be proportional to the energy spread as 

spread at the laser asymptotic equilibrium state. Jo IS the ft- t- t- bj t- 

,., .jT.; . 1. u 1 i-1 conventionally used in the accelerator community, 

spatial period between two electron bunches. Go the max- 

1 u li-.iiji .J These equations include the effect of the detuning by 

imum laser gam (which can be analytically determined ^ ° ■' 

i ^K ii ^a i jj J ii 1 i-.j-. r^i ji adding a cumulative delay between the electron bunches 
|251I28| and depends on the characteristics of the undula- ° •' 

, J 1 J- u ^ J r 4-1, .4-1 rpi and the laser pulse: 

tor and electron beam) and L the cavity losses, ihe syn- 

chrotron damping time Tg is related to the electron bunch t„+i = t„ -|- (5. (4) 
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This delay modifies the longitudinal overlap between the 
Gaussian electron bunch distribution and the laser pulse, 
which reduces the gain as: 



G„(r) = Go 



L 
Go 



cxp 



2al 



(5) 



The laser intensity at pass n can be expressed as the zero- 
order momentum of the laser intensity distribution : 

r'+oo 



yn{T)dT. 



(6) 



At perfect tuning, the laser tends towards an asymp- 
totic solution (constant pulse train shape), characterized 
by / = 1 and S = \ (see figure [T^), in agreement with 
model (|lallbp . When increasing the detuning, the laser 
intensity first remains constant in the stationary regime, 
then beyond a given threshold, it becomes pulsed (pulsed 
train shape) at a few hundred Hertz (see figure [T}d) . At 
this point, the resonant frequency Fjf changes with the 
detuning |29]. If one keeps increasing (5, a second thresh- 
old is crossed and the system quickly reaches a stationary 
state where both / and S are constant, and less than 
unity (see figure [T}:). It has been shown [16 that this 
model reproduces well the experimental dynamics of dif- 
ferent FEL facilities [5gi[5T|[52P5Pi] . For example, the 
absence of experimental stationary state at perfect tuning 
on the ELETTRA FEL has been explained by the small 
width of the central cw zone Note that similar dy- 
namics have been observed on conventional mode-locked 
lasers, which have been demonstrated to be governed by 
analogous equations [55] . 

In the following, the properties of the system are in- 
vestigated when the detuning is modulated at a frequency 




20 30 
Tirne (ms) 



4Q 50 



Fig. 1. Numerical simulations of the laser intensity versus time 
for different detuning values; a) 5 = fs, b) 5 = 6 fs, c) 5 = 36 
fs. Simulations parameters (corresponding to the Super- AGO 
FEL): TR = 120 ns, Go = 1.5 %, L = 1 %, ao = 6.5x10"'', 
(Te = 9.3x10-* , cjr = 107 ps, = 9 ms, is = 10"^ 

close to the resonant frequency of the system. Experimen- 
tally, the laser is first adjusted at perfect tuning (with a 
0.1 fs precision), then a modulation 5{t) is introduced with 
a function generator, which modulates the RF frequency 
at 100 MHs 



6{t) = asin(27r/t). 



(7) 



with / = 320 Hz. The control parameter a is tuned man- 
ually by adjusting the amplitude of the function generator 



^ Storage ring FELs are based on the emission of synchrotron 
radiation by charged particles in magnetic devices (such as 
dipoles or undulators), which consists in a conversion of ki- 
netic energy into a radiation. To compensate for the loss of 
energy at each interaction, and so at each pass in the cavity, 
the RF boosts the electrons: Its frequency has to be linked to 
the revolution period of the electrons in the ring. 
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signal. During the experiment, the laser intensity and its 
temporal spectrum (obtained by Fourier Transform os- 
cilloscope routine) are simultaneously recorded allowing 
changes in the regimes of the system to be detected. 

3 Characterization of the chaos 

In this section, the chaotic nature of the experimental FEL 
dynamics is explored, under a modulation of detuning. 
Analysis of the nonlincarities are performed on the ex- 
perimental signals. Figure [2] displays Super- AGO SRFEL 
macropulses, which have no apparent regularity and whose 
intensity peaks are strongly fluctuating. This is confirmed 
by inspecting the spectrum, which presents a continuous 
background between and 640 Hz (see figure [2Jd), which 
characterizes the aperiodicity of the time series. Nonethe- 
less, one can note that the spectrum exhibits a peak at 
twice the modulation frequency. 

In order to get a deeper insight into the dynamics, 
the resulting aperiodic time series are analyzed with the 
false nearest neighbours method to estimate the embed- 
ding dimension and with the calculation of the maximum 
Lyapunov exponent to appreciate the sensitivity to initial 
conditions [37] . 

Generically, introducing additional degrees of freedom 
into a system generates more complex dynamical regimes, 
since it can explore a phase space of larger dimension. 
Yet, systems with many degrees of freedom can exhibit 
macrostructurcs -both temporal or spatial- of a surpris- 
ing simplicity. In order to determine the effective com- 
plexity of a dynamics, one can resort, for example, to the 
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Fig. 2. Super- AGO chaotic dynamics for a detuning modula- 
tion of a = 16.5/s, at frequency 320 Hz. a) Laser intensity 
vs. time and its, b) Fourier Transform, with logarithm scale 
representation in the inset. The intensity measurements were 
performed with a photomultiplier (Hamamatsu R928), and an- 
alyzed by an oscilloscope (Lecroy LT264M, 350 MHz, l(f sam- 
ples/s) triggered by the trigger signal of the function generator 
(Stanford Research System DS345). The acquisition was set to 
25000 samples/s on a two-second time scale. The time series 
comprise at least 560 periods. 

embedding dimension, which gives an estimate of the num- 
ber of variables necessary to describe the dynamics. Such 
a reduction of the complexity is crucial, since in our case a 
continuous description of the light pulse distribution nec- 
essarily implies an infinite-dimensional phase space. The 
minimal dimension of the space in which a time series 



6 C. Bruni et al.: Chaos in 




□imaision 



Fig. 3. Ratio of false nearest neighbours versus the phase space 
dimension for three experimental time series of the laser inten- 
sity I{t) for a modulation frequency / = 320 Hz and a=33 
fs or 51 fs, calculated with the TISEAN programs [38ll39ll40j . 
using the delay coordinates. 

evolves can be estimated by the false nearest neighbours 
method :41j. One point I{ti) is "nearest" to another one 
/(ti+i) if the distance Ri between the two points is not 
higher than a heuristic threshold R^. Once the number of 
points for which Ri < Rt is zero or close to zero, then the 
embedded dimension is reached. 

Figure [3] relates the ratio of false nearest neighbours 
for three aperiodic time series of the laser intensity from 
Super-ACO experiment. This ratio falls to zero for a di- 
mension d = 3, which tends to indicate that the dynam- 
ics can be embedded into a three-dimensional space. De- 
spite of the complexity of the system, this result suggests 
that it may be possible to describe the dynamics in a low- 
dimensional model. 

The sensitivity to initial conditions can be character- 
ized by its Lyapunov exponent, that is, the divergence 
rate in phase space of two close trajectories [5711^ . Even 
though only the intensity can be measured, some tech- 
niques have been developed to rebuild a phase space from 
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Fig. 4. Calculation of In(^) as a function of time for different 
time series for a modulation frequency / = 320 Hz. The fit has 
been realized on the linear part. Programs of TISEAN [381139] 
has been used to draw this plot, using the delay coordinates. 

a single variable, such as the delay coordinates [/(i), I{t + 
(5t), I{t 2St),..., I{t + niSr)], where 5t corresponds to 
the delay. Although it provides a limited access to the 
dynamics, it is assumed that it can give a representative 
estimation of the chaotic properties of the system. 

Let us consider two close initial conditions distant by 
Aq in the rebuilt phase space. If we call At the distance 
between the states at ti + < and t2 + t, the Lyapunov 
exponent A/ is defined as: 

A; = lim lim log ^ (8) 

t^oozio^O Aq 

The Lyapunov exponent is then calculated in the phase 
space associated to the delay coordinates, with m = 2 in 
agreement with the result of the false nearest neighbours. 

Figure m represents the exponential growth in time of 
nearby trajectories. The linear increase relates the diver- 
gence, between 0.1 and 0.5 ms, where the linear fit al- 
lows the maximum Lyapunov exponent between 3.2 and 
3.8 ms~^ to be calculated. This positive exponent, esti- 
mated from experimental time series, is a quantification 
of the SIC of the Super-ACO FEL. 
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4 The route to chaos [S^3,b3,c3), and the intensity now oscillates at 2//3. Thus 

the laser response is now at //3 and its harmonics, al- 



In order to elucidate the emergence of chaos, previously 
revealed by the positive Lyapunov exponent for example, 
the modulation amplitude is progressively tuned. The ob- 
served dynamical regimes are analysed through their spec- 
trum and their trajectories in phase space. 

4.1 Dynamical regimes increasing the amplitude of the 
modulation 



though the even ones are at a much higher amplitude. The 
trajectory is a once again a closed curve, and its thickness 
in the (/ < l,dl/dt > 0) part attests to the presence of 
secondary peaks of intensity next to the main one. 

When the modulation amplitude is further increased 
until 35 /s, the laser dynamics becomes again chaotic. 
However one must note that although the trajectory fills 
an area of the reduced phase space similar to the one of 
the chaotic regime at 17 fs, the spectrum reveals that 
different frequencies are dominating the dynamics. 



Figure [5] shows five experimental dynamical behaviours 

representative of those observed during the different ex- 

; ,1 i • • • Finally, the modulation amplitude is set at 51 fs, and 

periments - there are some uncertainties m the experimen- ^ ' 

tal parameters. For a low detuning amphtude a = 7 /s, a ^^^t^"^ ^^^i^i*^ ^ ^^^^^^^ periodic regime, with the 

T-periodic regime is observed, caUed IT, with T = 1//. ^^"^ P^^^^^ intensity at a frequency, which is half the 

In this regime, the laser intensity (see figure EH) is mod- the a = 7 fs case. Once more, the presence of har- 

ulated at the modulation frequency (at / = 320 Hz), as ^^^^^^ frequencies can be associated with secondary peaks 

its spectrum shows. The projection in the (/, dl/dt) phase ^ ^^^'^^ intensity. 

space, called hereafter "reduced phase space" , is close to a Thus, the dynamics observed experimentally point out 

closed curve, which means that the system has reached a an alternation of chaotic and nearly periodic regimes, which 

hmit cycle. At a modulation amplitude of 17 fs, the laser enables the exploration of the route to chaos in the system, 

intensity exhibits an erratic behaviour: The delay between that is the transition from regular dynamics to chaotic 

the macropulses and the intensity peaks is now quite ir- ones. Nevertheless, a finer exploration of the bifurcation 

regular (see figure [5}2). Inspection of the Fourier space, diagrams remains experimentally problematic, because of 

where a broad spectrum is now observed, confirms this a drift in the system parameters. Actually, theoretical and 

aperiodicity. As the trajectory tends to fill a surface-like experimental studies on an internal modulated CO2 laser 

area in the reduced phase space, one can conclude on the have demonstrated that some bifurcations can disappear 

presence of chaos in the laser. because of perturbations present in the system [H]. 

Then, when the modulation amplitude is increased un- In order to further elucidate the mechanism of emer- 

til 28 fs, periodicity is back in the system (see figure gence of chaos in the system, the theoretical model is 
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Fig. 5. Experimental dynamics observed on the Super- AGO 
FEL. The laser intensity is displayed in the first column, its 
spectrum in the second one, and the projection of the trajec- 
tory in the reduced phase space in the last column. The so- 
called derivative coordinates [43 [I{t), dl{t)/dt, d? lit) /dt^ 
d"^I{t)/df^] are used to draw a representative picture of the 
phase space. The lines correspond respectively (from top to 
bottom) to the cases a = 7, 17, 28, 35 and 51 fa. Experi- 
mental conditions; / = 320 Hz, Ic = 30 mA, Go ~ 1.4 % and 
L ^ 0.8 %. 

used, since it allows for a finer tuning of the modulation 
amplitude. 

Because of the system parameter drift, a quantitative 
agreement between the experimental data and the model 
are hard to achieve. Instead, we choose typical values of 
parameters of the Super- AGO FEL (that is, 2% and 0.5% 
respectively for the gain and losses coefficients), and then 
we observe the succession of regimes which emerge, even 
if they occur at different modulation amplitudes. 

Figure [5] displays the regimes obtained from model 
I|3al3bp . At a low modulation amplitude (a = 24 fs), the 
intensity is regularly pulsed, at the frequency / (orbit IT). 
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Fig. 6. Dynamics of model pal3b|) . The first column shows the 
laser intensity time series, the second one its spectrum, and in 
the third one, the trajectory in the reduced phase space is 
represented. The lines depict, from top to bottom, the regimes 
at a = 24, 36, 108, 168 and 216 fs. The parameters are typical 
of the Super-ACO FEL: tr = 120 ns, Go = 2 %, L = 0.5 %, 
(To = 6.5x10"'', (Te — 7.5x10"* , (Jt — 107 ps, Ts — 9 ms, 
is = 10"^ / = 320 Hz. 

When a is increased to 36/s, the system switches to a 
chaotic regime, a peak at 2/ is still present, but the in- 
tensity spectrum is quite broad, and the trajectory begins 
to fill a mussel-like area of the reduced phase space. One 
has to note that the difference between the experimental 
and numerical trajectory in this space is essentially due 
to the measurement characteristic time, which is typically 
twenty times larger in the former case. Thus, a chaotic 
trajectory (such as case (c2) of Figs IS] and [H]), which is 
supposed to be dense in a given area of the reduced phase 
space appears to fill it very partially because of the finite- 
ness of the measurement time. Moreover, the shorter is 
this time, the smaller is the fraction of phase space that 
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it seems be filled. This behaviour is illustrated in FiglTl 
where two trajectories, measured over respectively 100ms 
and 2000tos, are compared. 
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Fig. 7. Experimental portion of the system trajectory in the 
reduced phase space, over a) 100 ms, b) 2000 ms. Experimental 
conditions: / = 320 Hz, a = 17 fs, h = 30 mA, Go = 1.4 %, 
L = 0.8 %. 

When a is increased from 36 to 108 fs, the intensity 
starts oscillating at frequency 2//3, and its harmonics (or- 
bit 3T): The system is now periodic, as it can be observed 
in the reduced phase space. The modulation amplitude is 
then set at a = 168 fs, a regime in which the intensity is 
composed of erratic peaks, the spectrum is broad, and on a 
time interval of 100 ms, the trajectory starts filling a sur- 
face of the reduced phase space, which is then a signature 
of chaos. Finally, when a reaches 216 fs, the intensity is 
pulsed periodically at //2, and the system seems to have 
a periodic behaviour. 

Thus, it was observed the succession of the IT-chaos- 
3T-chaos-2T regimes with the model (I3al3bp . which gives 
a qualitative agreement between the model and the experi- 
ments. Once more, the quantitative one remains problem- 
atic because of the system parameter drift in particular 
due to the electron bunches lifetime (which affects the 
number of electrons in the bunch and thus the gain) . Fig- 




a (fs) 



Fig. 8. Amplitude for which the transition of the dynamical 
regimes occurs as a function of the initial gain with fix losses. 
Same parameters as for figure (6) 

ure [8] shows that for different gain values, the succession 
of the regimes remains the same even if the amplitude for 
each transition changes with the gain. The model (|3a|3b|) 
reproduces well the alternation of periodic and chaotic 
regimes. 

4.2 Period-doubling bifurcation 

In order to elucidate the mechanism of emergence of chaos 
at a finer scale than the one allowed experimentally (not 
speaking of the parameters drift), the modulation ampli- 
tude is slowly increased in model pal3bp . using the pre- 
vious set of parameters. FiglTU] displays the evolution of 
intensity and of the system trajectory in the reduced phase 
space when a is increased. At a = 30 fs, the orbit is 2T. 
Indeed, although the intensity peaks still occur every 2//, 
their peak value is now modulated at half the modulation 
frequency. Then, an increase of 3 fs leads to a doubling 
of the period of the dynamics, that is an obit 4T: The 
intensity maxima are now modulated at //4. One more 
increase of 1.2 fs leads to an extra demultiplication of 
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Fig. 9. Simulated bifurcation diagram of the SRFEL: The 
intensity peaks are drawn versus the modulation amplitude. 
The cascade of periodic regimes eventually turns into chaos, 
before returning to periodicity. Same parameters as for figure 

El 

the orbit branches. And as seen previously, a little beyond 
that value (at a — 36 fs), the system dynamics becomes 
chaotic (see FiglS]). This period-doubling cascade is illus- 
trated in FiglQ] where the successive peaks of intensity are 
drawn versus the modulation amplitude. The doubling of 
the period of the regimes occurs at aiT-»2T ~ 27/s, then 
at a2T—>iT ~ 33/s, before the system enters a chaotic 
regime around a^x^c ~ 34/s. The return to a periodic 
2T regime is observed around ac~*2T ~ 44/s. New bi- 
furcations appear at higher values of a, which are not 
studied here. This study at constant gain/losses param- 
eters demonstrates that the Super-ACO FEL is entering 
a chaotic regime via a period-doubling cascade, a tradi- 
tional "route to chaos" (see fig.jS]). 




W 

limBims) l(ab.ijnit^ 

Fig. 10. Dynamics of model (|3al3bp for a fine increase of the 
modulation amplitude. In the left column is displayed the in- 
tensity evolution, in the right one, the trajectory in the reduced 
phase space. The three lines correspond, from top to bottom, 
to a = 30, 33 and 34.2 fs. Same parameters as for figure |6l 

For a higher ratio gain/losses (typically 3 times higher) 
than the experimental one, larger modulation amplitude 
step is required to go from orbit nT to 2nT. Since the 
step to go from 2T to 4T in the model is 3 fs, the ex- 
perimental observation of the doubling-period cascade is 
estimated around 1/s, which is less than the step applied 
for the modulation amplitude (1.2 fs). So that only the 
alternation of chaotic and regular regimes is observed ex- 
perimentally. The bifurcation diagram of the Super-ACO 
FEL can be compared to the bifurcation diagram of the 
logistic map. The first transition to chaos occurs through 
a period doubling cascade. Then other routes as intermit- 
tence should be observed, leading to periodic windows. 
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4.3 Intermittence 

A first manifestation of intermittence is the difference be- 
tween the two experimental chaotic behaviors at a = f 7/s 
and a = 35 fs (see figures]). The spectrum of the first 
chaotic motion has a dominant frequency at 2f/3 (see fig- 
ure [IHd). This is caUed hereafter 3C: C for chaos and 3 for 
3T. As a consequence, one can expect that a 3T-periodic 
orbit exists for close values of parameters, as it is effec- 
tively observed for a — 28 fs (see figurefTTbV By increasing 
the amplitude of the modulation, the spectrum blows up 
again, with a dominant frequency at f/3 (see figure fTTH ). 
This 3C evolves towards a chaotic regime with a dominant 
frequency at f/2, called 2C (see figure fTTb). The laser then 
becomes 2T periodic (see figure [TTf ) . A similar evolution 
of the attractor has been also observed with a modulated 
CO2 laser . This indicates that chaos and periodic dy- 
namics may occur through intermittence in a storage-ring 
FEL. A more detailed picture of the phenomenon can 
be found thanks to the model. In the simulated time sc- 
ries, burst of periodic signals appear in the chaotic signal 
(see fig. fT2|). The time between successive bursts tends to 
increase when getting closer to the periodic regime. The 
spectrum corresponding to a chaotic time series also ex- 
hibits a dominant frequency given by the periodic parts 
(see fig [T5|) . whose peak, at 2f/3, increases with the am- 
plitude of the modulation. When the intermittence is well 
established, the system switches to a periodic behavior 
whose frequency corresponds to the dominant peak in the 
adjacent chaotic windows. 
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Fig. 11. Spectrum of the laser intensity for different ampli- 
tudes of the modulation, a) IT behavior: a =7 fs, b) 3C be- 
havior: &— 17 fs, c) 3T behavior: a= 28 fs, d) 3C behavior: 
a= 35 fs, e) 2C behavior: a,— 50 fs, f) 2T behavior: a= 51 fs. 
Experimental conditions: Ic = 30 mA, Go ~ 1.4 %, L ~ 0.8 %, 
/ = 320 Hz. 

4.4 Symmetry of the system 

A manifestation of symmetry can be found in the time 
series given by the position of the center of mass of the 
light pulse called in the following as position. Figure [T3] 
shows the numerical results on the evolution of the laser 
pulse position. In all cases (periodic and chaotic) the evo- 
lution is periodically pulsed at the modulation frequency. 
The evolution of the amplitude of the position is well de- 
scribed by the projected phase space reconstruction from 
the derivation coordinates. For periodic cases, the phase 
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Fig. 12. Simulated laser intensity versus time with intermit- 
tence phenomena highlight in red. a=43.74fs. 

space shows an asymmetric limit cycle, whereas a sym- 
metric attractor is revealed for the chaotic regimes. In the 
system, there is a natural symmetry due to the fact that 
the modulation is symmetric compared to the zero tempo- 
ral detuning condition, leading to a pulse position which is 
naturally either positive or negative. These various sym- 
metries can be used to understand the chaotic dynamical 
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Fig. 13. Evolution of the intermittence with the modulation 
behavior, which opens new possibilities for a simplifica- amplitude and its consequence on the spectrum.l) Time series 
tion of the model. Note that the possibility for chaos of the laser intensity and its 2) power spectrum, a) a=44.04fs, 
exhibiting symmetries is stiU a debating topic and it has b) a=43.8 fs, c)a=43.74fs, d) a=43.68fs, e) a=43.62 fs. From e) 



been the subject of different studies [45ll46] 



5 Conclusion 



to a), the periodic outburst in the chaotic time series increases. 
Similarly, the spectrum shows a growth of the peak frequency 
at 2f/3 with f=320 Hz. 



In this paper, we have studied a particular case of elec- 
tron beam instabilities by applying an external modula- 
tion periodically modulated near the resonant frequency 
of the storage ring FEL system. It lead us to characterize 
the chaotic response of the Super-ACO FEL to such an 
electron beam instability. 

By calculating the embedding dimension of the time 
series, we demonstrate that a dimension three is sufficient 



to describe accurately the system phase space, including 
its nonlinear features, in spite of the infinite dimension of 
the model. The chaos is characterized by a positive Lya- 
punov exponent, which reveals the sensitivity of the time 
series to the initial conditions in the case of the Super- 
ACO FEL. The experimental data, combined with the 
finer analysis of the model, leads to the conclusion that 
the first transition to chaos is a period doubling cascade 
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Fig. 14. Numerical evolution of the center of mass of the FEL 
pulse for a modulation frequency f= 320 Hz. 1) Time series and 
2) reconstructed phase space of the moment of order two of the 
FEL pulse distribution, a) IT motion : &— 24 fs, b) Chaotic 
motion ; &— 168 fs, c) 2T motion : a= 216 fs. Numerical pa- 
rameters of figure IB] The FEL pulse center of mass is retrieved 
and analysed, through the first order momentum of the distri- 
bution. 

and the successive transition may occur by intcrmittence. 
The intermittence has been revealed experimentally by the 
evolution of the spectrum, where peaks could be observed 
as a sign of periodic orbits at close values of parameters. 
This signature of intermittence in the spectrum has been 
confirmed numerically. This analysis of the chaos reveals 
that a non conventional laser like the FEL has a bifurca- 
tion diagram very similar to the one of the logistic map. 
As a consequence, stabilisations of this kind of FEL in- 
stabilities can be envisaged and inspired for developing 
feedback systems, as it was already developed for table 
top lasers for example. 
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